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Abstract 



Numerical simulations of strongly compressible MHD corresponding to a stellar atmo- 
sphere with substantial gravity and near force-free magnetic fields show that the current 
sheet collapses (its width decreasing substantially). As a result, the reconnection rate in- 
creases dramatically 
> 
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I. INTRODUCTION 

It is well known that magnetic configurations are dissipating rather slowly in astrophysical 
conditions. Obviously, the Ohmic reconnection time, 

L 2 
rj 

where L is the characteristic magnetic scale, and rj is resistivity, is huge, comparable with 
cosmological times. 

The dissipation becomes much more efficient in the presence of current sheets. In par- 
ticular, the Sweet-Parker (SP) flow in the vicinity of the current sheet results in faster 
reconnection rate, so that it can be estimated as 

tsP = ^E = ^'\ (1) 



where S is Lundquist number, 



Ca is Alfven velocity, 



and tA is Alfven time, 
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tA = ± (2) 

see [I], |2f] . As S ^> 1 in astrophysical conditions, the SP mechanism (CQ) provides magnetic 
dissipation which is much faster than Ohmic. Still, this mechanism is rather slow compared 
with the time of stormy events like solar flares, or flares in the stars. 

More efficient mechanism was suggested by Petschek j3j], with reconnection time, 

81n^ 

tp = tA ■ 

IT 

The Alfven time (J2J) is usually small enough, so that the characteristic time of stormy and 
violent events is comparable or even slightly larger than it. For this reason, the time tp 
is sufficient to explain these stormy events. It was pointed out however, that there should 
be some additional nontrivial conditions satisfied in order this mechanism to work. In 



8) persistently showed that force-free 



particular, there are difficulties for this mechanism in simple two-dimensional geometry, 
without additional effects like anomalous resistivity, Hall currents, localized resistivity etc., 
see e.g. fl, fl, fl, |z|. 

II. IMPACT OF THE COMPRESSIBILITY. 

Most of the numerical simulations and analysis are provided by incompressible or only 
slightly compressible situations, although it was suggested long ago that compressibility may 
speed up the reconnection rate, [9|. It was indeed shown experimentally that the compression 
plays significant role in the reconnection process, [10 [ 

Previous quite extensive numerical simulations 
compressible reconnection rate increases dramatically. These simulations also confirm that 
compressibility plays key role in the reconnection process. In particular, it was noticed that 
the current sheet is strongly compressed, or one may say that it collapses, increasing the den- 
sity. As a result of compression, the sheet width shrinks, leading to increased reconnection 
rate. 

The compression is substantial in the low density plasmas, typical for force-free (or near 
force-free) solar and stellar atmospheres, where the forces are balanced not by a regular 
pressure, but rather by magnetic pressure. 

The plasma and current sheet collapse is in fact forming a singularity, suggested by B.C. 
Low, |Uj]. That was strictly one-dimensional model, and in fact, in the vicinity of the current 
sheet the magnetic configuration is indeed quasi-one- dimensional. 

The simplest way to see the plasma compression is to consider an one-dimensional model, 
where all the quantities depend on x only, and magnetic field B = {0, B y (x,t),B z (x,i}. 
The force-free equilibrium can be sustained if 

B y + B l = COIlst - ( 3 ) 

We suggest that B y is an odd function of x, changing sign at x — 0, then B z is an even 
function, and (if B z > 0) with maximum at the origin. Suppose that the characteristic 
length of By is 5, then, according to ([3]), the characteristic scale of B z would be the same. 



Assuming that 5 is "small" (we will define this scale later), we conclude that B z has a sharp 
maximum at the origin. 

The equation of B z evolution is as follows, 

d t B z + V • vB z = V V 2 B Z (4) 

The velocity consists in fact of only one component, v x , and, as seen from (jlj), it is an odd 
function vanishing at the origin. Consider then (jlj) at x = 0. If we look for stationary 
solutions, d t = 0, and, as v x (x = 0) = 0, we get, 

B z (x = 0)V • v = r]d x d x B z (5) 

As the right-hand-side is negative (corresponding to the maximum), we conclude that V- v < 
0, i.e., compression. We now suppose that 5 correspond to the width of the SP current sheet, 

S = 5sp = -±- r (6) 

Then, for not very strong B z , corresponding to the compressible media, we have, 

*—*"£■ (7) 

It follows from (JTj), that, in case of SP current sheet, the compression proceeds with the 
fastest time in the problem. 

The meaning of this compression is as follows. As mentioned above, the equilibrium 
([3]) can be sustained if magnetic pressure B^/(8n) has a sharp maximum at the origin 
compensating magnetic "wall" consisting of sharp increasing of By. Without this maximum 
of B z the configuration would collapse to the x = point (due to the By/(8ir) pressure). 
However, finite diffusivity smoothies out this maximum, not allowing the equilibrium ([3]) to 
survive. This plasma collapse tries to involve magnetic field due to the frozen-in conditions, 
thus resulting in magnetic collapse as well. It is this situation that was observed in previous 
simulations and also will be reported below. 

III. THE ROLE OF THE GRAVITY 

Theoretically, we can consider force-free fields (of stellar coronas), completely neglecting 



the pressure, as in Uj. In that case, the compressible part of the Lorenz forces are balanced 



by magnetic pressure. In fact, any compression in plasma should increase the magnetic 
pressure gradient to balance and eventually to stop the compression. However, previous 
numerical simulations persistently showed that magnetic pressure is not sufficient, and that 
the regular pressure grows. As follows from Sec. [TTl the magnetic pressure B^/(8n) is 
unable to balance all the forces because of the finite diffusivity, quite substantial in SP 
current sheets. 

In spite of the growing regular pressure, we did not seem to find anything dramatic in 
its impact. The reconnection is still much faster than in SP current sheets. In particular, it 
follows from ([1]) that the reconnection rate for SP flows is 

r S p = ~ V lf \ (8) 

tsp 

and, as rj is low for astrophysical conditions (in dimensionless units), the reconnection is in- 
efficient. In previous simulations, with pressure included, the reconnection rate r ~ ^o.2h-o.3^ 
i.e., much faster than the SP rate. 

Still, the pressure seems to somewhat slow down the reconnection. We note, that dramatic 
increase of density perturbation, and therefore increase of the pressure, is unlikely in the 
stellar atmospheres. Indeed, unless there are special conditions (satisfied, e.g., for the solar 
prominences), in the presence of gravity forces, the density "blobs" cannot be sustained in 
coronas for a long time, and they are supposed to slip down along the magnetic field lines. 

We expect that, if the fall down time is much less than the Alfven time, then, roughly 
speaking, the matter remains in stratified and almost unperturbed state all the time during 
the reconnection process. 

We introduce a dimensionless number, 

gL gL 



where g gravity acceleration, 



CW2 " C A s : 



a p 



(9) 



and s is the sound speed. Note that we consider low pressure case, i.e., /3<1. 

In equilibrium state, the matter is stratified along the field lines. We suppose that the 
pressure is constant on the photosphere, and therefore, in gravitational equilibrium, this 



stratification is independent of field line position, i.e., the density and pressure are functions 
only of vertical dimension. In other words, in equilibrium state, the forces 

V x B x B 

gp-v P + — (10) 

are balanced separately, i.e., there is hydrostatic equilibrium of a stratified media, 

g p-Vp = 0, (11) 

and magnetic equilibrium, i.e., force- free field, 

V x B x B = 0. 

These two equilibrium's have two different time-scales to establish. We will look for the 
case when the first equilibrium has a shorter time-scale, and therefore established, while the 
magnetic field is evolving slower. 

The estimated scale-height of the stratified media (ITTj) is 

ff 4 (i2) 

We will suppose that 

Z > (3 l /\ (13) 

a condition, easily satisfied in the stellar atmospheres (like the sun). Then, 

H < L. (14) 

Then, the recovery time for any density perturbation to vanish and to return to the strati- 
fication fTTTT) would be just the falling time on one scale- height (TT2|) . 



m V2.9 
t/ = VT = — ' (15) 

Obviously, we are interested in the case of fast established stratification, that is, 

t f < t A , (16) 

which is equivalent to the requirement 

Z > 1, (17) 
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an inequality stronger than (1131) . 

In this paper we are dealing only with the limiting case (IT7I) . leaving more general case 
for the future studies. Presenting density and pressure as 

p = po + 5p, 

P = Po + 5p. 
Thus, in the zeroth approximation, the first two terms in (fTOl) cancel each other, 

gPo - Vp = 0, 

cf. (ITT]) . Assuming that (fT6l) . or f|T7|) is satisfied, the density perturbation vanishes for a 
short time tf, being generated for a much longer time t^, according to ([7]). Therefore the 
density perturbation is weak, 

In the first approximation, the gravitational force is no longer balanced by the gradient of the 
pressure, although they are of the same order of magnitude. Both terms can be estimated 
as 

r PO 



while the estimation of the Lorentz force is 

Po 



r<2 



L 
Then, the ratio of the first two terms in (1101) to the Lorentz force is 

-^- = — = B^ « 1 

Summarizing: In the zeroth approximation, the first two terms in ( llOp are canceling each 
other, and in the first approximation they are small compared to the Lorentz force, if 
Z>1. We will consider the momentum equation in this approximation. As the buildup of 
the density and therefore of the pressure, according to (TTHj) . is weak, we do not expect the 
pressure to impede the reconnection process, as opposed to our previous simulations. 



IV. FORMULATION OF THE PROBLEM 

The problem is solved in two-dimensional (compressible) MHD. That means that all 
quantities depend on x, y only (and independent of z), but the fields contain all three 
components (the so-called 2.5-dimensional MHD). Thus, 

B = {B ± (x,y,t), B z (x,y,t)} = {d y A z (x,y,t), - d x A z (x,y,t), B z (x,y,t)}. (19) 

The initial configuration is depicted in Fig. [TJ The magnetic field is rooted on the photo- 
sphere, y — 0, decreasing dramatically up into the corona. The density is stratified, i.e., 
it is decreasing (exponentially) from the photosphere into the corona. The configuration is 
considered to be periodic in horizontal direction (i.e., in X-direction). We are solving the 
induction equation, 

9,B = Vx[vxB]+ r?V 2 B, (20) 

with zeroth initial velocity v. The momentum equation is presented by Navier-Stokes equa- 
tion with forcing as in fflOl) . in approximation of Z 3> 1, see Sec. IIHI 

We can imagine that two "sunspots" of opposite polarities are initially far away of each 
other, and therefore they are not connected. Suppose that, in the course of evolution, 
they approach each other, forming a configuration similar to that depicted in Fig. [TJ Due 
to frozen-in conditions, the topology of the field lines is not easy to change, and, for a 
while, they remain disconnected from each other. Due to finite conductivity, a current 
sheet will form, and the reconnection starts. Further evolution is observed in our numerical 
simulations. 

The case of B z = corresponds to current free equilibrium. That means that, if the 
pressure is neglected (which is justified in low-/? plasma (of stellar atmospheres), the only 
equilibrium to which the configuration is tending to settle down is current-free. Obviously, 
the configuration depicted in Fig. (TJ does not correspond to this equilibrium, and it has 



to form a current sheet, the situation described by Syrovatsky, 12|. If B z ^ 0, smooth 
distributed currents (not current sheets) may be formed in corona, so that the configuration 
may settle down to force-free equilibrium with (smooth) currents in the corona. This situ- 
ation is much less understood, being actually typical for stellar atmospheres, and therefore 
is interesting for numerical studies. 
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FIG. 1: Initial configuration 

The Alfven velocity, is small on the photosphere, and it reaches a maximum in the middle 
corona, and again very small (numerically zero) is in the upper corona. For this reason, 
almost nothing happens on the photospheric level, where the Alfven velocity is 100 times 
less than its maximum, and, for the same reason, practically everything is still well above 
in corona. The reconnection is expected to take place where the Alfven speed is maximal, 
and it is indeed observed in the vicinity of its maximum. 

Naturally, we have in mind a real three-dimensional configuration that can be approxi- 
mated as a two-dimensional. That means that the field in ignorable direction z is eventually 
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also closing at the photosphere. In self-consistent approach, we may assume that the nor- 
mal component of the velocity vanishes on the photosphere, which means, in particular, that 
v z = 0. 

The configuration depicted in Fig. [TJ or other similar configurations with slightly different 
geometry, were treated in several (at least three) different codes, with quantitatively similar 
results, described in the tho next sections. 

V. RESULTS OF NUMERICAL SIMULATIONS: MEASUREMENTS OF THE 
FIELD EVOLUTION 

Previous studies were devoted to the so-called rosette-structure, in which case the current 



sheet formation is inevitable, see 



It has been shown that the reconnection rate strongly 



increases, as compared with the SP rate, 8|, see also the main results in Fig. [HJ All previous 
simulations are dealing with no gravity, Z = 0, see definition in (Q. The new simulations are 
devoted to the opposite limiting case, Z ^> 1, and they are performed with pseudo-spectral 
code. 

A. Dramatic magnetic field reconstruction 

The main evolution of the magnetic field is described by the behavior of the vector- 
potential, defined in ffT9l) . According to (120]) . the governing equation for it reads, 

d t A z + vVA 2 = V V 2 A Z , (21) 

which is equivalent to the ^-component of the Ohm's law, 

E + vxB = ??j, (22) 

where 

d t A z = -E z . (23) 

In astrophysical conditions, dimensionless value of r\ <C 1, and therefore the right-hand- 
side in (122]) . is negligible, that is, in fact, the electric field is defined from 

E + v x B = 0, (24) 
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FIG. 2: The reconnection flux versus time; r/ = 0.00125, for fully nonlinear case from numerical 
simulations (solid line), and for pure Ohmic decay (dashed line). 

which corresponds to the frozen-in conditions. This is true everywhere except at the x-point 
where B x = B y = 0. Here, the typical electric field is weak, and, as seen from (1221) . is defined 
only by finite (and very weak) Ohmic decay. 

The situation is dramatically different in the presence of strong currents, - current sheets, 
which appear in the vicinity of the x-points. It is in these places where the reconnection 
occurs. In our case, the x-point occurs somewhere on the V-axis, possibly moving up and 
down along the axis in the middle corona. 

In initial configuration, Fig. (TJ A z (x = 0, y) = on the F-axis. Afterwards the reconnec- 
tion occurs, and A z grows on the axis. Almost nothing occurs both on the photosphere, at 
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y — 0, and high in the corona, where the Alfven velocity is small. As a result, A z (x = 0,y) 
acquires a maximum corresponding to the x-point (where B x = B y = 0). The value of the 
maximum, A™ ax , characterizes the reconnected flux. It is important to estimate, how fast 
is this maximum growing, that is, how much flux is reconnected during a fixed amount of 
time. In fact, the time derivative of A z in this point, according to ff23l) . defines the electric 
field (responsible for the particle acceleration to high energies, etc.). In some runs, the max- 
imum of electric field may present the efficiency of the reconnection process. However, this 
happens only if this maximum is not sharp, i.e., it is surrounded by high values of electric 
field - in the vicinity of the maximum. Otherwise, the maximum is short-lived. That means 
that relatively fast reconnection did happen, but only for a short time, so that only a small 
fraction of magnetic flux reconnected. For this reason, we will calculate average electric 
field, which is compared with maximum possible electric field, 

£max = Q A B U (25) 

see (El]), where the velocity was replaced by Alfven velocity, maximal in this problem. 

Figure [2] illustrates one of the runs of our simulations. It can be seen that, for only several 
Alfven times, a substantial amount of magnetic flux is reconnected. This is compared with 
pure Ohmic decay, when all nonlinear effects are switched off. Note that Ohmic decay 
(with superconductive boundary condition at the bottom and zero at the top, and periodic 
conditions on the left and right), also results in reconnection, but very slow, as seen from 
the figure. 

The final state of magnetic field is depicted in Fig. [HI which should be compared with 
initial, Fig. [TJ It can be seen that much of the flux has been reconnected. Note that, 
if all the nonlinear effects are switched off, that is, only Ohmic decay remains, then, the 
configuration at the same moment of time would be practically the same as initial depicted 
in Fig. [TJ 

The magnetic energy is decreasing quite efficiently during the evolution, see Fig. HJ It is 
clear that the corresponding energy drop during the same time in linear case, that is, only 
due to pure Ohmic decay, is much less efficient. The energy at t = 5, corresponding to the 
configuration in Fig. [3], proved to be close to the theoretical final energy, as seen from the 
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FIG. 3: Final configuration, after reconnection; t = 5 and r/ = 0.00125 



Fig. H] (dashed-dotted line). 



B. Strong currents 



Initial magnetic field is not at equilibrium, and, of course, the system is trying to reach 
equilibrium state. This can be achieved either by current-free configuration, for small or 
absent B z , or force- free state, - for finite B z . However, the initial topology of the field lines 
is incompatible with final equilibrium state, see 12j |. and therefore the field lines should 
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FIG. 4: Evolution of magnetic energy, for rj = 0.00125, Z ~^> 1, depicted with a solid line. This is 
compared with pure Ohmic decay (dashed line) , which would eventually end up in the lower state 
depicted with dashed-dotted line. 

reconnect. As the system is driven to equilibrium by Lorentz forces, the reconnection should 
be faster than Ohmic. This can be achieved if strong currents are generated, i.e., in fact 
current sheets. 

As seen from Fig. strong currents are indeed formed. An equilibrium state would be 
presented by smooth currents, like the initial in the figure, but this state is achieved only 
after part of the flux is reconnected, and this process of reconnection is accompanied by 
strong currents. 
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FIG. 5: Section of currents at x-point at t = 1, r\ = 0.00125, solid line with asterisks. It is compared 
with initial current - smooth dashed line. 

C. Strong compression 



As mentioned above, the equilibrium state can be achieved only after some flux is recon- 
nected. Thus, the current sheet is inevitably formed. Of course, the starting width of the 
sheet is defined by SP mechanism, see (E]). In the vicinity of the sheet the configuration is 
in quasi-one-dimensional equilibrium ([3]). As mentioned in Sec. [TTJ i? 2 -component acquires 
maximum at the center of the sheet. For this, there should be some compression, that is, 
5 z -component should increase approximately in 
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FIG. 6: Density buildup in case of rosette structure. Initial density is unity 

times, - in this rough estimate we use initial characteristic values. This means that the 
density would grow the same way. For previous simulations of rosette- structure, this estimate 
was 1.2-r 1.5, or so, while the real compression proved to be much stronger, see Fig. [61 The 
mechanism describing further compression after "natural" due to the maximum of the B z - 
component was described in Sec. [Til In addition, it can be seen from Fig. [6] that the 
compression grows with decreasing resistivity with some scaling explained in Sec. IVIAI 

For Z ^> 1, the real evolution of the density is defined by the gravity forces, and, as 
mentioned in Sec. IIH| is not significant. However, we can calculate the "pseudo-density", 
that is, knowing the velocity field from the simulations, we can calculate the density which 
would be real if the gravity force are switched off. Using Lagrangian solution of the mass 
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FIG. 7: Density buildup at the Y-axis, depicted by solid lines. Dashed lines correspond to the 
initial state. For both cases rj = 0.00125, and (a) is for t = 1 and B±/B z = 100, while (b) 
corresponds to t = 5, and B±/B z = 2 

conservation equation, we obtain results some examples of which are depicted in Fig. [3 

It is obvious that, for strongly compressible case, depicted in Fig. [7(a), the compression 
is much stronger than "natural", that is only 100, while in essentially incompressible case, 
depicted in Fig. [7J(b), the compression is not significant. We note however, that the be- 
havior of "pseudo-density" is only an illustration of strong compression, as opposed to the 
measurements of the density evolution in previous simulations with Z = 0. 
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VI. RESULTS OF NUMERICAL SIMULATIONS: THE SCALING 

A. Ohmic decay and scaling relationships 

The process of reconnection is accompanied by Ohmic decay. In regular astrophysical 
conditions it is quite weak, because the field is frozen-in into the plasma. The situation is 
different in the presence of current sheets. The energy release during the Ohmic decay is 
described as follows, 



0_ 
Ft 

or, as an estimate, 



I B 2 dxdy = -rjM = -rj /"(V x B) 2 dxdy, (27) 



r B " L2 « ^5L, (28) 

tR d 2 



where on the right-hand-side the volume of the current sheet of the width 5 is estimated as 
SL because, from the numerical experiment, the length of the current sheet does not change 
much and it is taken to be L. 

Then, the estimated reconnection rate, 



1 rj 

WR = 

or, in dimensionless form, 



UR -T R -W 



For SP current sheet, 5 is defined in (jBJ), and, according to (1251) . 

Vsp = ^, (30) 

cf. (0Q). Further compression of the sheet width, that is, decreasing of 5, would result in 

increase of the reconnection rate. In order to estimate this compression, we write ( l29l) in 

equivalent form, 

J^6sp_ l_ _ Ssp 1 
R ~ 5 SP 5 S~ 5 S 1 / 2 ' 

or, the compression can be defined as 

£ = 5 ~f = ^^ 1/2 , (31) 
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FIG. 8: The reconnection rate, exponents with power-law fitting. Previous runs, Z = 0, /3 = 0.01, 
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and 2>1 (diamonds) is essentially incompressible. 



where po is initial density. This relationship between the reconnection rate and compression 
will be used below. 

Note that direct measurements of the current sheet width 5 could be bias. Much more 
reliable are measurements of integral quantity M defined in (J2"T|) . the measure of Ohmic 
decay. According to 

(32) 



B 2 I 2 

M « Q R S 
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FIG. 9: Ohmic dissipation, M, denned in (|27p . as a function of Lundquist number S. Here Z 3> 1, 
and B z /B± = 0.01. 

B. Scaling laws 

Our main results are summarized in Fig. [SJ Previous calculations correspond to Z — 0, 
and they definitely show reconnection faster than SP rate. All new simulations correspond 
to the opposite limiting case Z > 1. If B z /B± = 0.4 then the velocity proved to be 
incompressible, and then the SP mechanism is recovered (depicted with diamonds and fitted 
with short dashed line). On the other hand, if B z /B± = 0.01, the plasma is strongly 
compressible, as illustrated in Fig. [3 Then, the scaling (fitted with a solid line) is practically 
flat. Note that the asterisks in the figure correspond to the measurements of averaged electric 
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field, see Subsection IV A\ normalized on maximal electric field, f l25|) . Therefore, the actual 
numbers should be compared with unity. It can be seen from the Fig. [S]that the numbers 
are reasonable, not very small. 

Another check of the consistency of the theory is the scaling of the density build-up 
depicted in Fig. [61 According to (13T|) . the compression should scale as 5' ' 23±0 ' 05 , which is 
consistent with the observed scaling given in Fig. [6j 

Finally, the scaling of Ohmic dissipation is shown in Fig. [91 This is roughly consistent 
with predicted in ( J32i) . although it does not exactly correspond to it, presumably because 
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all observed scaling laws can be trusted up to the first decimal number, in spite of the fact 
that the dispersion of the points is less than that, as can be seen in all figures. 

C. Scaling laws as a function of compressibility 

We also provided numerical simulation in the same scaling range as in Fig. [8] for differ- 
ent compressible media, while all other initial parameters (like density and magnetic field 
distribution, etc.) are the same as in all simulations. The compressibility is changing with 
the parameter (B z /B±) 2 . When this parameter is small, then the media is highly compress- 
ible, while increasing this parameter, we will approach essentially incompressible situation. 
As seen from Fig. (TUJ the incompressible media is reach when (B z /B^) 2 = .04, and then 
the Sweet-Parker approximation is recovered. For higher compressibility we found that the 
reconnection is faster (i.e., the scaling exponents are smaller). Note that previously simula- 
tions correspond to moderate compressibility, in which case the reconnection is faster than 
SP mechanism, but still slower than what we have in case of high compressibility. 
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